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NELREPORTIISS 


Study  vlf  and  If  eleGtromagnetiG  propagation  meGhan-^ 
isms.  Study  the  interrelationships  of  phenomena  within 
the  Upper  stratosphere  and  ionosphere  whiGh  affect  vlf  and 
If  eleetromagnetie  propagation  and  the  extent  to  which 
variations  of  these  phenoinena  may  influenee  radio  propaga^ 
tion.  The  study  reported  here  examined  the  correlation 
between  the  minimum  virtual  J’^layer  height  h ' E  at 
Washingtoni  Di  C.  ^  and  the  stratospheric  circulation  over 
the  eastern  United  States  as  a  function  of  frequency  and 
scale  size  of  the  circulation. 


RESULTS 

The  minimum  virtual  J*- layer  height  was  found  to  be 
correlated  with  a  particular  type  of  large’-scale  circulation 
at  the  lOO^mb  level  which  has  a  period  of  6,  3  days.  One  of 
the  ten  circulation  indices  used  to  desGribe  the  100*^mb 
circulation  was  found  to  have  predominant  periods  of  22 
and  12  days.  The  predominant  12 -day  period  appears  to 
be  controlling  the  correlation  between  the  100-mb  circular 
tion  and  h,' E  at  the  6, 3^day  period. 


RECOMMENDATION 


Correlate  other  ionospheric  parameters  with  the 
circulation  over  larger  regions  of  the  world  and  at  other 
pressure  levels  in  the  stratosphere  or  troposphere. 
Examine  the  physics  of  the  relationship  between  h'^  and 


the  stfatosphefic  circulation  particularly  with  respect  to 
the  preference  for  h*E  to  be  related  with  certain  circular 
tion  patterns  which  occur  at  certain  periods. 
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INTRODUCTION 

Many  studies  have  been  direeted  toward  the  Gommon 
goal  of  deterniining  the  relationship  between  the  tenuous 
ionosphere  and  the  considerably  denser  atmosphere  below, 
the  stratosphere  and  tropospherei  (See  list  of  references 
at  end  of  report. )  This  study  examines  the  relationship 
between  the  minimum  virtual  J’=layer  height  h ' E  and  the 
lOO^mb  circulation,  as  a  function  of  frequency  and  scale 
size  of  the  circulation. 


DEFINITION  OF  VARIABLES 


h'E  is  a  function  of  the  temperature  structure  of  the 
atmosphere,  solar  radiation,  atmospheric  composition, 
and  recombination  rates.  h'E  represents  the  lowest  region 
of  the  ionosphere  regularly  measured  and  is  essentially  the 
height  of  the  base  of  the  I'^layer,  usually  at  about  100  km, 
h'E  at  Washington,  D.  C. ,  was  chosen  to  represent  the 
ionosphere  and  was  taken  from  the  CRPL-F  series  pub¬ 
lished  by  the  National  Bureau  of  Standards,  ^  The  height 
is  tabulated  in  kilometers  for  every  hour  during  the  day 
when  the  layer  exists.  The  error  of  b'Ev^as  ±10  km 
before  January  1956  and  ±2  km  afterwards.  These  large 
errors  mask  most  variations  thought  to  be  of  interest,  but 
they  were  effectively  reduced  by  taking  a  seven-hour  aver¬ 
age  h'  oih'E  centered  around  local  noon  each  day,  h ‘ 
should  be  a  better  daily  value  of  h'E  than  a  single  hourly 
value. 


It  was  decided  to  examine  the  relationship  between  the 
circulation  on  a  constant  pressure  surface  and  h'  because 
both  are  functions  of  the  atmospheric  temperature  struc¬ 
ture.  The  lOO-mb  winds  at  the  seven  stations  shown  in 
figure  I  were  chosen  to  represent  the  lOO-mb  circulation. 


Ftguti  2*  The  •#v«n  ttgtipnt  U99d  to  d99erth9  tho 
lOO^mb  eirouldtion* 


The  number  of  stations  was  limited  by  the  eapacity  of  the 
eomputer  for  the  statistical  techniques  chosen  a.nd  the 
sample  size  of  2922  days.  The  station  separation  was 
seleeted  SO  that  relatively  larger  and  small'seale  sizes 
of  the  Gireulation  Would  be  included.  The  10O*mb  winds 
were  supplied  by  the  National  Weather  Records  Center  for 
the  eight^year  period  from  March  1951  through  February 
1959, 


The  wind  directions  were  given  to  16  points  of  the 
compass,  and  the  speeds  were  given  in  meters  per  second. 
The  daily  and  u*  coniponents  of  the  winds  were  computed 
for  each  station.  The  u-'component  is  positive  toward  the 
east  and  the  y-component  is  positive  to  the  north. 
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The  u^-  and  u“ components  and  h'  were  filtered  to 
remove  the  annual  eyele  in  the  meanSi  The  filter  response 
given  in  figure  2  shows  that  the  amplitude  reduction  of  the 
annual  cycle  was  0*  05  while  the  amplitude  of  all  variations 
with  periods  less  then  100  days  was  unaffected.  The 
sample  size  for  the  filtered  data  was  2722  days^ 


rtgurt  Betponae  eturva  of 
tho  mathomat ieal  ftltor  uaod 
to  ftltor  tho  data» 


It  is  possible  that  h'  is  related  to  eertain  scale  sizes 
of  the  rOd»inb  circulation  at  certain  frequenciesi  A  factor 
analysis  technique  developed  in  the  Department  of  Mete¬ 
orology  at  the  Massachusetts  Institute  of  Technolo^^  was 
used  to  generate  circulation  indices  which  represent 
various  scale  sizes*  The  technique  produces  what  has 
been  labeled  "empirical  orthogonal  functions.  "  The  tech¬ 
nique  was  used  to  generate  14  imcorrelated  time -varying 
weighting  factors  3^(1  stsl4)  of  14  corresponding  orthog¬ 
onal  circulation  patterns  y-i  from  the  correlated  daily 
values  of  the  seven  pairs  of  filtered  u-  and  u-  components. 
This  technique  factorized  the  total  variance  of  the  correlated 
u-  and  u-  components  into  a  minimum  number  of  orthogonal 
variables  3^'s  which  contain  a  certain  per  cent  of  the  total 
information.  The  3^'s  and  y. 's  make  up  the  "empirical 
orthogonal  functions.  "  Those  g^'s  containing  little  variance 
can  generally  be  considered  to  be  noise  with  respect  to 
application  in  prediction  and  need  not  be  correlated  with 
h'.  Out  of  the  14  2^'s  generated  from  the  14  u-  and  u- 
components,  ten  g^'s  having  the  largest  variance  were 
required  to  contain  90  per  cent  of  the  information.  The 
remaining  four  g^'s  with  small  variance  were  eliminated 
in  the  correlation  analysis.  The  first  ten  g^'s  were  the 
circulation  indices  used  in  this  study.  The  distribution  of 
the  variance  among  the  14  g^’s  is  shown  in  figure  3, 


Figura  3»  Tha  diatrihut ion  of  tha  oari^ 
anea  among  tha  id  g^^a* 


The  U“  and  u*  components,  the  2j's,and  the  y^'s  are 
related  according  to  the  matrix  equation: 


U-  V-  -  .  i  .V 
17  1117 


Tlf  Tl\  717 


/^ll  “ ‘^i  .lol  /u^ii  “‘U^17  ~vhi  •• 


^  *ni  ^n,i( 


•  •  • 


u^i  o  i  1  *  ‘  ‘U  ^1 6 , 7  y^i  6  ,  i 


where  n  is  2722  days.  Each  row  in  the  U‘v  matrix  is  the 
set  of  and  u*  components  for  the  seven  stations  on  day  t 
Each  column  i  in  the  2  ^matrix  is  a  function  of  time  t  and 
is  orthogonal  to  the  other  nine  columns.  Each  row  t  in 
the  y-matrix  is  a  function  of  space  and  is  orthogonal  to  the 
other  nine  rows.  The  u-  and  u-^components  for  station  J 
On  day  t  are 

is 

^tj  ~  \  u^tj 

t-1 


ip 

i“i 

From  equation  2  the  direction  and  magnitude  of  all  winds 
can  be  determined  with  the  noise  removed,  K  all  S^'s  are 
zero  except  when  i  ^  s,  then 


iH=  kts{u^sj*ysj)^^ 

tans  .  ^ 

^  uVsJ 

When  f  g  is  negative,  the  wind  direction  is  opposed  to  the 
direction  given  by  equation  3b, 


Figure  4  shows  the  first  ten  orthogonal  funGtions  of 
spaGe  derived  from  equations  3a  and  3h.  The  isopleths 
are  relative  height  contours  in  feet  for  equal  to  1,  It  is 
noticed  that,  in  general,  the  scale  size  of  the  eireulation 
deGreases  as  i  increases.  This  agrees  with  earlier  studies. 
It  is  now  evident  that  a  correlation  between  h'  and  the  g^'s# 
the  weighting  factors  of  the  y^'s,  is  a  correlation  between 
h‘  and  various  scale  sizes  of  the  lOO^^mb  circulation. 


Y,  Y 

Fiffurt  Th*  ftrmt  t«n  orthogonal  funotiona  of 

apaeo,  Tho  iaoplotho  ara  rolativo  hotght  contour# 
of  tho  lOO'^mb  ourfaoo  whon  q.  ^  1»  Tho  wind 
opoodo  givon  at  oaeh  mtation^ aro  in  motora  par 
aoeond* 
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POWER  SPECTRUM 


If  the  relationship  between  h '  and  the  g  ^  Vs  is  to  be 
examined  as  a  funGtion  of  frequency  or  period,  as  well  as 
a  function  of  scale  size,  cross  ^spectrum  analysis  between 
and  the  g^'s  is  required^  Cross  “spectrum  analysis 
gives  the  square  of  the  correlation,  the  coherence  Oti,  as 
a  function  of  period,  and  the  phase  lag  t  between  correlated 
frequencies.  Two  time  series  can  be  related  at  certain 
frequencies  which  are  not  predominant  in  either  time 
series.  An  examination  of  the  power  spectrum  oth'  and 
the  g^'s  will  indicate  the  relative  amplitudes  of  all 
frequencies  present  in  the  time  series. 

Figure  5  shows  the  uimormalized  power  spectra  ufyi) 
of  h'  and  the  first  ten  g^'s  as  a  function  of  the  ordinate 
number  n  for  n  between  0  and  80.  The  frequency  interval 
between  ordinates  is  equal  to  4  x  lor^  cycles  per  day. 

Thus  the  periods  range  from  infinity  to  3,  3  days.  There 
are  43  degrees  of  freedom  in  each  estimate  of  u(n).  All 
power  spectra  generally  show  a  concentration  of  energy  at 
low  frequencies  and  a  gradual  decrease  of  energy  as  fre^ 
quency  increases,  a  common  feature  of  meteorological 
data.  The  concentration  of  variance  in  the  lower  order 
g^'s  should  be  noted. 

No  predominant  periodicities  are  present  in  the  power 
spectra  except  the  spectrum  for  g^ ,  There  are  two  very 
significant  periodicities  in  gg  as  indicated  by  the  90  per 
cent  confidence  limits.  *  The  two  predominant  periods  are 
12.  5  and  22. 2  days.  A  power  spectrum  of  g^  with  20 
degrees  of  freedom,  which  is  not  shown,  mdicates  an  even 
greater  level  of  significance  at  the  two  periods  and  places 
the  predominant  periods  at  12,  5  and  21,7  days.  The 
presence  of  two  predominant  periodicities  in  g^  implies 


*The  true  value  of  Uij%)  will  lie  within  the  limits  given  in 
figure  5  for  90  per  cent  of  the  time  when  there  are  43 
degrees  of  freedom. 


that  |/g  varies  every  12*  5  and  22  days  from  a  ey clonic  to 
an  anticyclonic  circulation  as  gg  goes  from  a  negative  to 
a  positive  value*  The  amplitudes  of  the  two  predominant 
periods  are,  of  course*  additive*  No  pressure  pattern 
motion  is  implied  because  pg  is  a  fixed  pattern  of  space* 

If  yg  is  part  of  a  wave  pattern  around  the  Northern  Hemi  = 
sphere,  then  pg  is  part  of  a  standing  wave  pattern*  y, 
extends  about  3000  km  in  the  east^west  direction.  If  pg  is 
one-half  a  600d*km  wave,  then  pg  suggests  a  wave  number 
of  5*  Based  on  a  numerical  experiment,  Phillips^  pre* 
dieted  that  the  wave  number  5  would  have  an  eastward 
velocity  of  1800  km/day*  It  is  interesting  to  note  that  the 
beat  between  a  12*  S  and  21,7  day  period  is  a  29*  5  day 
period,  which  is  almost  a  lunar  cycle* 


Figure  The  power  mpgetra  of  h*  and  the  firat 

ton  orthogonal  funotiona  of  timt  with  43  degreoa 
of  frtodow,  Tha  fraauanoy  interval  betwaon  tha 
ordinataa  ia  4  x  io~^  cyoiea/day*  Tha  amplitudaa 
of  tha  low  fraguano iaa  hava  baan  raduoad  by 
filtering  the  wind  eomponanta  faaa  fig*  2)* 


Ftgurw  5«  (Oont tnutd) 


Utn» 
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CORRELATION  BETWEEN  It' 

AND  THE  lOO  mb  dROILATION 


Figure  6  shows  the  coherenGe  OH  between  h'  and  the 
first  ten  g^’s  for  periods  between  250  and  5  daySi  The  OH 
between  two  series  of  eight-digit  random  numbers  having 
the  same  degrees  of  freedom  is  shown  for  eomparisom 
Although  independenGe  between  the  g^'s  at  all  periods  has 
not  been  established,  the  OH's  of  all  ten  g^'s  were  added  to 
give  the  total  OH^  as  shown  at  the  bottom  as  0H^»  If  the 
ten  g^'s  are  orthogonal  at  all  periods*  then  OH^  is  equiva¬ 
lent  to  the  total  reduction  in  varianGe*  oH^  can  be  used  as 
a  guide  in  determining  the  most  important  periods  until  the 
degree  of  independenGe  between  the  g^'s  has  been  estab¬ 
lished* 

There  are  five  predominant  peaks  in  the  ten  OH  funGtions 
in  figure  6.  Two  peaks  oceur  with  g^*  one  ooeurs  with  gg, 
and  two  oGGur  with  q^,  g^  is  related  to  h'  at  a  period  of 
about  13  or  14  days  and  at  a  period  of  6.  3  days.  The  shape 
ofthe  peak  regions  in  OH  for  g^  and  g^  are  surprisingly 
similar  between  the  6.  3  and  6. 6  day  period,  g^  also  has 
a  peak  in  OH  at  a  period  of  17.  9  days.  The  probability  that 
the  five  peak  values  of  OH  could  happen  by  chance  ranges 
from  0. 165  to  0.  035  as  shown  in  figure  6.  Although  the 
five  OH  peaks  are  not  highly  significant,  an  examination  of 
the  relationship  between  h',  g^,  g^,  and  g^  at  a  period  of 
6,  3  days  was  thought  to  be  valuable  primarily  because 
the  relationship  apparently  involves  the  fundamental  periodic 
ity  of  12,5  days  found  to  be  predorninant  in  g^  .  OH^^  also 
indicates  that  the  relationship  between  h '  and  the  1 00-mb 
level  is  singificant  at  a  period  of  6.  3  days.  The  remaining 
part  of  this  paper  will  be  an  examination  of  the  relationship 
between  h',  ,  q^  ,  and  g^  at  a  period  of  6,  3  days. 

The  phase  lag  t  in  days  for  q^,  q^ ,  and  gg,  relative  to 
h'  for  a  period  of  6.  3  days,  is  2,  4,  4,4,  and  6.2  days 
respectively.  The  90  per  cent  confidence  limits  of  these 


PERIOD  T  (DAYS) 

125  62.5  41.7  31.5  25  20.8  17.8  15.6  13,9  12.5  11.4  IG.4  9:6  8.8  8.3  7.7  7.4  7.0  8.6  6.5  6.0  5.7  6.4  5.2  50 
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figure  6,  Th«  eohertnee  OH  hftwttn  the  firet  ten 
9^'#  and  h*ae  a  funetion  of  the  ordinate  humher  n 
or  period*  The  moot  a ignifioont  regiono  of  OB 
hove  been  aooentod* 
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lags  are  *0^  7,  40*  8*  and  40,  7  days  respeGtively.  *  None  of 
the  90  per  cent  limits  overlap,  can  be  assumed  to  be  in 
phase  With  h'  because  the  lag  of  6.  2  days  is  essentially 
equal  to  the  period  of  6.  3  days.  Assuming  that  the  phase 
lags  are  2*  4,  and  0  days  for  ^  *  and  g^  ,  respectively, 

and  that  the  relative  amplitudes  at  a  period  of  6,3  days  are 
given  by  the  power  spectra,  the  cyclical  pressure  pattern 
related  to  h '  can  be  constructed. 

Figure  7  gives  the  cyclic  pressure  pattern  related  to 
h'  as  specified  by  g^,  gg,  and  gg  for  an  even  six*day  period 
instead  of  a  6.3 “day  period.  The  isopleths  in  figure  7  are 
relative  heights  of  the  100~mb  surface  for  the  relative  wind 
speeds  shown.  On  Day  1  (4)  a  trough  (ridge)  extends  in  the 
east“west  direction.  On  Day  2  (5)  the  trough  (ridge)  begins 
to  tilt  in  a  NE“SW  direction.  On  Day  3  (6)  the  southern 
part  of  the  trough  (ridge)  moves  eastward,  while  the  norths 
ern  position  remains  about  stationary,  but  the  general  NE* 
SW  tilt  is  conserved.  The  maximum  (minimum)  in  h ' 
Occurs  on  Day  3  (6),  The  pressure  trough  on  Day  3,  when 
K'  peaks,  is  quite  similar  to  a  wave  on  a  frontal  surface. 
Apparently  is  associated  with  NE*SW  oriented  ridges 
and  troughs  on  the  lOO-mb  surface,  which  occur  about 
every  six  days.  The  maximum  of  h'  occurs  when  gg  is 
negative  and  increasing;  that  is,  when  the  cyclonic  flow  is 
decreasing  and  being  replaced  by  an  anticyclonic  flow.  The 
phase  lag  between h,'  and  g^  at  a  period  of  about  13  days  is 
about  10  days,  which  indicates  that  the  maximum  Of  h' 
Occurs  when  g^  is  about  zero  and  decreasing  in  magnitude. 


The  true  value  of  t  will  lie  within  the  stated  limits  90 
per  cent  of  the  time  when  there  are  43  degrees  of  freedom. 


PAY  3  PAY  4 


figure  7,  The  eix-’-dag  <iyci/c  IQO^mb  preeeure 
pattern  that  ie  related  to  h**  Tho  ieoplethe  are 
relativo  height  contour#  of  tho  lOQ'^ub  eurfaoe. 


Figure  ?•  (Continued) 


100“  mb  Q‘ CIRCULATION 
DECEMBER  29,1956 


lOO  “  m  b  0  “  CIRCULATION 
DECEMBER  26,1956 


100  -  mb  Q  “  CIRCULATION  100  mb  Q  •  CIRCULATION 

DECEMBER  30,1956  DECEMBER  31, 1956 


Figure  9.  The  100-^mh  preaaure 
pattern  from  28  Deo  ember  1956 
to  1  January  1957  accord ing  to 
the  firat  ten  q^'a  given  in 
figure  8*  The  iaopletha  are 
relativa  height  contoura  for 
the  lOO^mb  aurfaca. 


IQO-mb  Q. CIRCULATION 
JANUARY  1, 1957 


COHERENCE  BETWEEN  THE  CIRCUUTiON  INDICES 

To  determine  the  Use  of  CH^  as  an  indicator  Of  the  total 
reduction  in  variance,  (7^  between  andgg*  and^g, 
and  fg  and  were  determined.  These  QH's  are  shown  in 
figure  10*  It  should  be  noted  that  relatively  large  peaks 
Occur  in  CH  at  a  period  of  about  6.  3  days,  particularly 
between  32  and  gg  and  between  jg  andgg*  It  should  be 
noted  that  neither  gg  *  nor  gg  has  a  predominant  peri^ 
odicity  at  a  period  of  6.  3  days*  The  peak  between  gg  and 
gg  has  a  p -value  of  0*  052  and  is  more  significant  than  the 
other  peaks*  The  similarity  in  the  shape  of  the  qH  curves 
of  h'  ‘■q^  and  h'-gg  at  a  period  of  about  six  days  probably 
resulted  from  the  correlation  between  gg  and  gg  at  about 
a  six-day  period.  Thus,  it  appears  that  OH^  can  be  used 
only  to  indicate  the  most  important  regions  where  % '  and 
the  lOO-mb  circulation  are  related. 


Figure  10,  Th»  eoh0r9na9  CB  bftvffn  9^  and  qg» 
Qg  and  qg,  and  qg  and  qg* 
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CONCLUSIONS 


From  the  standpoint  of  statistiGal  significaneei  the 
most  important  result  of  this  study  is  the  deteetion  of  two 
predominant  periodicities  in  one  of  the  ten  orthogonal 
indices  used  to  describe  the  lOO^mb  eirGulation.  These 
predominant  petiods  occur  in  and  are  22  and  12,  5  days. 
Apparently  h'  and  the  lOO^mb  circulation  are  only  related 
at  periods  of  about  18, 13,  and  6  days*  The  correlation 
between  h'  and  at  an  IB^day  period  was  not  examined  in 
this  Study.  The  meehanisms  producing  the  22  and  the  12.  5 
day  period  in  g^  niay  be  different  because  h'  is  not  related 
to  the  I00=mb  circulation  at  the  22*^day  period  while  it  is 
at  the  12.  5  “day  period.  It  is  thought  that  the  statistically 
deduced  pressure  pattern  related  to  h'  at  a  period  of  about 
6*  3  days  is  significant  although  the  significance  of  the  OH 
peaks  at  a  period  of  6.  3  days  is  not  high  because  the  6.  3“day 
pattern  is  apparently  controlled  by  the  predominant  12. 5 “day 
period. 


Correlate  other  ionospherie  parameters  with  the  eireula= 
tion  over  larger  regions  of  the  world  and  at  other  pressure 
levels  in  the  stratosphere  or  troposphere.  Examine  the 
physics  of  the  relationship  between  h'E  and  the  stratospheric 
circulation  particularly  with  respect  to  the  preference  for 
h  'E  to  be  related  with  certain  circulation  patterns  which 
occur  at  certain  periods. 
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